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We consider non-Abelian dyons in Einstein- Yang-Mills-Higgs theory. The dyons are spherically 
symmetric with unit magnetic charge. For large values of the electric charge the dyons approach 
limiting solutions, related to the Penney solutions of Einstein-Maxwell-scalar theory. 
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I. INTRODUCTION 



Magnetic monopoles arise as regular non-perturbative solutions of SU(2) Yang-Mills-Higgs theory Monopoles 
with magnetic charge n = 1 are spherically symmetric, monopoles with higher charge have axial symmetry 3, [H, HI, HI 
or no rotational symmetry at all fGj. When electric charge is added to magnetic monopoles, dyons arise 0, Q. These 
Y^' are stationary solutions with vanishing angular momentum [H, [l(| ■ 

The nontrivial vacuum structure of ST/ (2) Yang-Mills-Higgs theory not only allows for magnetic monopoles, but it 
^Q 1 results in a plethora of further regular non-perturbative solutions [ll|, [l2|> [l3|, [l5j . The simplest of these solutions 
are monopole-antimonopole pairs, where a monopole and an antimonopole form an unstable equilibrium configuration 

When gravity is coupled to Yang-Mills-Higgs theory, gravitating monopoles [la . [T7| . dyons [18|], monopole- 
antimonopole pairs [H, H3, HH, [H[ , and further configurations arise [23[ ■ These solutions depend on a dimensionless 
coupling constant a, which is proportional to the square root of Newton's constant and the Higgs vacuum expectation 
, value. 



For each type of solution, a branch of gravitating solutions emerges smoothly from the corresponding flat space 
solution and extends up to a maximal value of a, beyond which the size of the core of the solution would be smaller 
than its Schwarzschild radius [lj| . At the maximal value of a this fundamental branch bifurcates with a second branch 
of solutions. 

For monopoles and dyons, in the case of vanishing Higgs potential, this second branch reaches a critical value of a, 
where it bifurcates with the branch of extremal Reissner-Nordstrom black holes with the same charge(s) @, Oil [13] • 
For monopole-antimonopole pair solutions and other composite solutions, in contrast, the second branch extends back 
to a — 0, where a pure Einstein- Yang-Mills solution is reached (after scaling w.r.t. a) [24. l25l. l26l |27|. 

As the second monopole resp. dyon branch reaches the critical value of the coupling constant a, and merges with 
the extremal Reissner-Nordstrom branch, the spacetime splits into two regions [161 ]. The exterior spacetime of the 
critical solution, corresponds to the one of an extremal Reissner-Nordstrom black hole [lQ], while the interior spacetime 
retains regularity at the center, due to the influence of the non-Abelian fields present. 

Here we reconsider gravitating dyons and study the properties to the solutions in the limit of large electric charge 
Q. We demonstrate, that for Q — > oo we also obtain a limiting solution, which consists of two regions, a non-Abelian 
interior regions and an Abelian exterior region. But in this case, the non-Abelian interior part has a flat metric, and 
the Abelian exterior part corresponds to a certain Penney solution of Einstein-Maxwell-scalar theory (after scaling 
w.r.t. Q) H|. 

In section II we present the Einstein- Yang-Mills- Higgs action, the Ansatz, the equations of motion and the boundary 
conditions. We then consider the equations for large electric charge, and relate them to the equations of Einstein- 
Maxwell-scalar theory. We describe the relevant features of the Penney solutions in section III, and discuss the 
properties of the numerically constructed non-Abelian dyons in section IV, where we also address their relation to the 
Penney solutions. We present our conclusions in section V. 
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II. EINSTEIN- YANG-MILLS-HIGGS SOLUTIONS 
A. Action 



We consider the SU(2) Einstein- Yang-Mills-Higgs action in the limit of vanishing Higgs potential, 



S = 



-gd A x (l) 



16ttG 2 v ^ '4' 

with curvature scalar i?, SU(2) field strength tensor 

F/u/ = dyA, - duA^ + ie [Ap,A v ] , (2) 

gauge potential A^ = l/2T a ^4°, gauge covariant derivative 

D M =V M + ie[A ll ,-] , (3) 

and Higgs field <& = r°<I> a ; G is Newton's constant, and e is the gauge coupling constant. Since we here consider 
vanishing Higgs potential, we impose a Higgs field vacuum expectation value v. 

Variation of the action Eq. fl} with respect to the metric g llUl the gauge potential A" and the Higgs field <& a leads 
to the Einstein equations and the matter field equations, respectively. 

B. Ansatze 

We employ Schwarzschild-like coordinates and parametrize the line element by 

ds 2 = -A{rfN{r)dt 2 + -^dr 2 +r 2 dn\ M r )=l-^, (4) 

N(r) r 

with mass function fi(r). The Ansatz for the gauge potential and Higgs field is given by 



A^dx* 1 = - — (r v dQ - t 6 sin 0dip) + ^p- T rdt , $ = u$i (r)r r , ( •". ) 



where the su(2) matrices i>, re, and are defined as scalar products of the spatial unit vectors with the Pauli 
matrices r° = (T x ,T y ,T z ), and the subscripts on the functions indicate the correspondence to the functions of the 
more general Ansatz, necessary for the construction of monopole-antimonopole systems [bl. [25j. 

C. Equations of Motion 

We now change to dimensionless quantities, the dimensionless coordinate r, the dimensionless electric function B\, 
and the dimensionless mass function p,, 

r = evr , B\ = — - , ft = evfj, , (6) 
ev 

and introduce the dimensionless coupling constant a 

a 2 = AttGv 2 . (7) 

Suppressing the ~ in the following, the tt and rr components of the Einstein equations yield for the metric functions 
the equations, 

M' = « 2 (^f + + + \Nr 2 ^ 2 + a^l! + ^ , (8) 

and 

2 / 2B 2 H 2 2K' 2 , 2 \ , 
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while the matter field equations yield 



(ANH' 2 )' = AH 2 [B-l + - JL ] , (10) 

^y = « (11) 

and 

(r 2 A/V$i)' = 2A$i# 2 . (12) 

These equations then depend only on the dimensionless coupling constant a 116]. 

A particular solution of these equations is the embedded Reissner-Nordstrom solution with unit magnetic charge, 
P = 1, and arbitrary electric charge Q, 

p(r)=m- a2(Q 2 a r +1) , A(r) = l, (13) 

ffa(r) = 0, Si(r)=i/-^, *i(r) = 1 , (14) 

r 

where the extremal solution satisfies 

m = r H = a^Q 2 + 1 • (15) 



D. Boundary Conditions 

Dyons are globally regular particle-like solutions of the SU(2) Einstein- Yang-Mills- Higgs system. Regularity of the 
solutions at the origin then requires the boundary conditions 

/i(0) = , (16) 

and 0] 

ff 2 (0) = l, Bi(0) = 0, *i(0) = 0. (17) 

Asymptotic flatness of the solutions, on the other hand, implies that the metric functions A and fx both approach 
constants at infinity. We adopt 

A(oo) = 1 . (18) 

The matter functions satisfy the asymptotic boundary conditions 

# 2 (oo) = 0, Bi(oo) = v, $i(oo) = l. (19) 

The dimensionless magnetic charge P is the topological charge of the solutions, while the dimensionless electric 
charge Q, and the dimensionless scalar charge ch are obtained from the asymptotic expansion of the fields, thus 

P = 1 , Q = — lim r (P>i — v) , c H = lim r 2 d r <S>i , (20) 
and ^(cxd) = m represents the dimensionless mass of the solutions. 
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E. Scaled equations of motion 

Since we are interested in the limit of large electric charge, Q — > oo, we now consider the above set of equations in 
terms of quantities scaled by the charge Q. Thus we introduce 

r m _ c H ( . 

r= Q' m= Q' CH= Q- (21) 



*' = ?W^ + tf*?) . ( 22 ) 



In the limit Q — > oo we then obtain a coupled set of equations for fi, A, B\ and $1, 

2 \A 2 

A' = a 2 f$?A , (23) 

L r = Q = 1 > (24) 

f 2 AiV$i = c H , (25) 
where the ' denotes differentiation w.r.t. f, while the equation for H2 decouples, 

H 2 " = H2(H ; 2 - 1] . (26) 

The coupled set of equations Eqs. P2"|) - (f2l)|) corresponds precisely to the set of equations of Einstein-Maxwell-scalar 
theory, when the electric charge has the value Q = 1. We therefore now turn to the Einstein-Maxwell-Scalar Solutions. 

III. EINSTEIN-MAXWELL-SCALAR SOLUTIONS 
A. Action 

The matter Lagrangian of Einstein-Maxwell-Scalar theory reads 



£ = -jF^F^- ia^iS"*! (27) 



with the Abelian gauge potential A^dx* 1 and the scalar field $1, giving rise to the stress-energy tensor in the Einstein 
equations, 



Gfj.v = 2a 2 T^ , (28) 



in terms of dimensionless coordinates and fields. 



B. Penney solutions 

Penney obtained static spherically symmetric solutions with electric charge Q and scalar charge ch- He employed 
the metric parametrization 



ds z = -e~ a dt 2 + e a dR 2 + e b dn 2 , (29) 



with radial coordinate R. 

The Penney solutions then read 



.■" = !!/?-/;ii/?--"ii '( ^ ^ — ) * >'■ = (SB-, MR -«))<" . (SO) 

a — p 
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e b B[ = Q , e 6 - a $i = c H , (31) 

where Aq = B\ denotes the time component of the gauge potential. 
The constants p, a, and A of the Penney solutions satisfy the relations 

A > = a 2 Q 2 , o?cl = (1 - A 2 ) (E^j 2 . (32) 

Since we want to relate the Penney solutions to the limiting dyon solutions, we must select those Penney solutions, 
which do not exhibit a metric singularity anywhere except the origin. Thus we now focus on the case, where in 
the Rcissncr-Nordstrom limit no horizons occur, but a naked singularity is present. In the Penney solutions this is 
achieved by the choice p = po + iao, a = po — iao- Defining 

f = (R - p ) 2 + o% , tan*=-^— , tan^o = -— , 7 2 =Po + ^o, (33) 

it — po (To 



then yields the metric functions 



cos(A* - tpo) 

COS (po 
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e b = £ 2 e a (34) 



and the relations, Eq. (|32l) . 

a 2 Q 2 a 2 c 2 
A 2 = — f- , cos 2 (po = 9 ^, 9 H 9 , m 2 = a 2 Q 2 sin 2 (p , (35) 

where the mass m is obtained from the asymptotics. From sin 2 tpo + cos 2 (po = 1 we obtain 

1 r 2 

= 1 ■ (36) 

Q 2 -(?) 

Since A 2 > we obtain the bound 

4<Q 2 -(^) 2 - (37) 



C. Limit of large electric charge 

We now consider the limit of large charge, taking Q — > oo. Then according to Eq. (|3"5|) also A — > oo. Eq. (|3^|) then 
requires that the bound in Eq. (|37[) is precisely saturated, 

4 = Q 2 -(^) 2 - (38) 

To obtain the limiting Penney solution for Q — > oo, we thus impose this bound. We then consider the set of 
equations w.r.t. the scaled coordinate R = R/Q, and introduce the scaled quantities Q = Q/Q = 1, 7 = 7/Q, and 
ch = crIQ- In particular, we reexpress the relation a 2 Q 2 = j 2 A 2 as a 2 Q 2 = 7 2 A 2 , i.e., 7^0. In this limit 




We are interested in the outer extremum of the metric function e a , since this is a particular point of the spacetime. 
It will later be identified as the transition point, where the spacetime of the dyons will split into a non-Abelian interior 
region and an Abelian exterior region. 

The outer extremum of e a occurs at -Ro, where 

^S=arctanji--1, e a ^ = ± . (40) 
ito y c H c H 
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We now impose the dyon boundary conditions on the Maxwell and scalar function, 

B 1 (R )=0, B 1 (oo) = u, *iORo)=a, $1(00) = 1. (41) 
Integration of the scalar field Eq. (f3"Tj) then yields 

poo 1 - 

^(oo) = c H / — dfi=-jE = l, (42) 



i.e., -Ro = c H , which together with Eq. (|4*D|) leads to 

cos a = c H , (43) 

relating the scaled scalar charge to the coupling constant a. 

To integrate the gauge field Eq. (|3ip . we first reexpress the metric function e° via 



cos 


a ( 


C H 








R 





cos a 



(44) 



Integration then yields a relation between the asymptotic value of the gauge potential and the coupling constant a, 

- , , f°° h - sin a 

Bi(oo) = / e~ b dR= = 1/. (45) 



i.'u 



a 

Addressing finally the transformation to Schwarzschild-like coordinates, we note that 

A 2 N = e- a , r 2 = e b , dR = Adr . (46) 
The transition point ro is thus given by r$ = e b ( Ro ' = RqQ 2 /^ = Q 2 , i.e., fo = 1. 

IV. RESULTS 
A. Gravitating dyons 



Gravitating dyons have been considered before [13, |22| • We here reconsider gravitating dyons and address their 
dependence on the electric charge Q. We focus on large values of the charge and, in particular, the limit Q — > 00. 

For dyon solutions in flat space, the magnitude of the electric charge is uniquely determined by the asymptotic 
value v of the electric component of the gauge potential. This potential parameter v is bounded, < v < 1, and in 
the limit of vanishing Higgs potential, one has the monotonic relation between Q and u, 

Q(v) = ^JL=, (47) 
V 1 — v 

Thus the charge diverges in the limit v — > 1. 

For gravitating dyons it was also expected, that (in the limit of vanishing Higgs potential) the bound v = 1 would 
be reached monotonically and that it would correspond to infinite electric charge |la |. 

When constructing the dyon solutions numerically, however, a surprise is encountered. For fixed coupling constant a, 
the gravitating dyons do not vary monotonically with v. Instead at a maximal value u max a bifurcation is encountered. 
Here a second branch extends slightly backwards towards smaller values of v. These bifurcating branches are exhibited 
in Fig. [TJ where we demonstrate the dependence of the charge Q and the scaled mass rh = m/Q on the parameter v. 

As seen in Fig. Q] the second branches rise very steeply as Q becomes large. They are confined to the intervals 
fcr(oi) < v < ^ m ax(«)- The lower bounds z^ cr (a), approached in the limit of infinite charge, are very close to the upper 
bounds v mEUt (a). The critical values v c -[(a) are exhibited in Fig.[2j We note, that for small a, the critical values p C i{<x) 
exhibit an almost quadratic dependence on a. No solutions are found beyond a = tt/2. 

Let us now consider dyons for very large values of Q, in order to identify the solution obtained in the limit Q — > 00. 
Clearly, the mass m diverges in the limit Q — * 00. However, the scaled mass fh = m/Q tends to a finite limiting 
value, which depends on the coupling strength a, as seen in Fig. [TJ Another quantity of interest is the scalar charge 
ch, determining the 1/r power law decay of the Higgs field of the dyon. Like the scaled mass, the scaled scalar charge 
ch = ch/Q approaches a finite limiting value, depending on the coupling strength a. In Fig. [5]ch is shown versus a 
for a very large value of the electric charge, Q = 10000. 



7 



(a) 



(a) 



100 
80 
60 
40 
20 



80 



40 



a=0.5 



0.959 0.9595 



a=0.0 
a=0.3 
a=0.5 
a=0.8 
a=1.0 
a=1.2 



a 



0.5 



0.6 



0.7 



0.8 



0.9 



0.241 
0.24 



0.9582 



a=0.5 



0.9592 



a=0.1 
a=0.3 
a=0.5 
a=0.8 
a=1.0 
a=1.2 
a=it/2 




0.2 



0.4 



0.6 



0.8 



(b) 



FIG. 1: Dyons: (a) electric charge Q(v) versus v, (b) scaled mass m = m/Q versus v, for several values of a. 
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FIG. 2: Dyons: (a) critical values v cr versus a, (b) scaled scalar charge ch = ch/Q versus a. The exact dependence as obtained 
from the Penney solution is also exhibited for v C r and ch. 

B. Relation with the Penney solutions 



To understand the limiting behaviour, we now consider the solutions themselves for very large values of Q. We 
exhibit in Fig.[3]thc metric and matter functions in scaled Schwarzschild-like coordinates, f — r/Q, for a dyon solution 
with very large electric charge, Q — 10000, at a coupling strength a = 0.5. The solution then appears to consist of 
two parts, an interior part in the region < f < fo, and an exterior part in the region fo<f< oo, with fo = 1. 

In the interior region, < r < 1, the limiting solution is given by 



A(f)= const, N(f) = 1 , Bi(f) = $i(f) = , H 2 = H 2 (f) 



(48) 



i.e., the metric is flat and the time-component of the gauge potential and the Higgs field both vanish. The only 
non-trivial function is the spatial gauge potential function H 2 , satisfying the single decoupled equation, Eq. (|26[) . 

In the exterior region 1 < r < oo, on the other hand, the spatial gauge potential function H 2 vanishes, while the 
metric and the other two matter functions satisfy the coupled set of equations Eqs. (|22p - (l231) . This system represents 
a special case of the coupled Einstein-Maxwell-scalar equations, where the scaled electric charge has the value Q = 1, 
studied in a different parametrization by Penney (28[. 

As discussed in section III, in the relevant Penney solutions no horizons occur, but a naked singularity is present. 
In the limit Q — > oo, these Penney solutions precisely saturate the bound on the scalar charge, Eq. (|57|) . i.e., they 
satisfy a quadratic relation between the scaled global charges, 

to 2 = a 2 (1 - cfj) . (49) 

Due to the dyon boundary conditions for the Maxwell and scalar functions the Penney solutions furthermore satisfy 
the relations Eq. (|4"5|) and Eq. (|45[) . i.e., the scaled scalar charge ch, and the potential parameter v are given in terms 
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(a) r (b) r 

FIG. 3: Dyons: (a) metric functions A and N, (b) gauge potential and Higgs field functions Bi, H2, $1, versus scaled 
Schwarzschild-like coordinate f — r/Q, for an almost limiting solution at a — 0.5, Q = 10000. Also shown is the limiting 
solution: the Penney solution in the exterior, and the flat non-Abelian solution in the interior. 



of the coupling strength a, 



sin a 

ch = cos a , v = , (50) 



and thus the scaled mass satisfies fh = a sin a. The validity of these relations for the limiting dyon solutions is seen 
in Fig. [T]for the scaled mass fh = m/Q and in Fig. [5] for v and ch = cjj/Q. 

Identifying the transition point fo as the outer extremum of the metric functions gu and g rr of this particular 
Penney solution, we find fo = 1 (Eq. PS|) ). Evaluating the metric functions A and iV at the transition point fo = 1 
yields 

A(f ) = ch , iV(fo) = 1 . (51) 

For comparison we superimpose in Fig.[3]the limiting solutions in the two regions. In the exterior region f < f < 00. 
the limiting solution is the Penney solution with the same a and Q = 1, saturating the mass bound, Eq. (|49|) . This 
Penney solution also determines the constant metric functions in the interior region < f < fo, via their boundary 
values at the transition point fo, Eq. (|5ip . On the other hand, the boundary conditions at the origin determine the 
constant matter functions B\ and $1 in the interior, and thus also at the transition point fo- 

In contrast the matter function H2 is a solution of Eq. (J^U) [29] . Expanding the solution H2 (f ) of the Yang-Mills 
equation Eq. (121H) in a power series, i?2(f) = 1 + h%f+ hif 2 /2 + 0(f 3 ), shows that h\ = and hi is a free parameter, 
characterizing the solution. Solving the equation numerically in the interval < f < 1 with boundary conditions 
#2(0) = 1, ^2(1) = h* , and varying h* , we observe that Ii2(h*) increases with increasing h* and tends to a finite 
value /i' 2 nax ~ 3.047 when h* tends to infinity. Comparison with the numerical dyon solutions for large charge Q, we 
find some evidence that H 2 (0) indeed tends to the value /i™ ax for Q — > 00, i.e., H2 diverges at f = 1 in the limit. 
Convergence is very slow, however, approximately like 1/VQ- 

We finally consider the range of a where dyon solutions exist. For small electric charge it was shown before [Tc| , 
that dyons exist only below a maximal value of a, e.g., a max (<3 = 0) = 1.40 and ce max (Q = 1) = 1.41. As the charge 
is increased further, a max increases as well. In the limit Q — > 00, we obtain from the above considerations a bound 
for a, 

< a < I . (52) 

Concerning the limit a — > tt/2 we see that the scaled scalar charge ch tends to zero. Thus in the limit there is no 
(non-trivial) scalar field in the exterior, whereas there is still an electric field. This suggests, that for a — > ir/2 the 
scaled limiting solution corresponds in the exterior to an extremal Reissner-Nordstrom solution with charge Q = a. 
Inspection of the analytical formulae and the numerical solutions shows, that this is indeed the case. We exhibit the 
approach towards this limit for the metric and matter functions in Fig. [4] We note, however, that in this limit, the 
RN solution is approached only in the interval 7r/2 < f < 00. The limiting solution is discontinuous at 7r/2, since 
the metric function A(f) assumes the RN value A = 1 only for f > 7r/2, whereas it vanishes for f < 7r/2, causing the 
metric to differ from the RN metric in this region. 
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FIG. 4: Dyons: (a) metric functions A and N, (b) gauge potential and Higgs field functions Bi, H2, $1, versus scaled 
Schwarzschild-like coordinate f — r/Q, for solutions with Q — 10000 and several values of a, approaching a — n/2. Also shown 
is the limiting solution: the extremal RN solution in the exterior, and the flat non-Abelian solution in the interior. 



Concluding, we see that the scaled limiting solution consists of two parts: an Abelian exterior Penney solution, 
saturating the bound Eq. (|49p and determined by the value of a (with Q — 1), and a non-Abelian interior part with 
flat metric, where the constant metric functions are determined by the exterior Penney solution at the transition point 
f = l. 



V. CONCLUSIONS 



We have reconsidered dyon solutions of Einstein- Yang- Mills- Higgs theory, in the limit of vanishing Higgs potential. 
Dyon solutions then exist for arbitrarily large values of the charge. As the charge becomes large, the coupled system 
of Einstein- Yang-Mills-Higgs equations decomposes into a single equation for the magnetic gauge potential, and a 
coupled system of equations for the metric, the electric gauge potential and the scalar field. This system of equations 
is equivalent to the set of Einstein-Maxwell-scalar equations, studied by Penney in a different parametrization [28l ] . 

In the limit Q — ► 00, the non-Abelian dyons solutions then tend to limiting solutions, which consist of two parts, a 
non-Abelian interior solution with flat metric, and a gravitating Abelian exterior part. The interior solutions represent 
non-trivial solutions of the single equation for the magnetic gauge potential (with the remaining equations trivially 
satisfied). They are regular at the origin, but diverge at the transition point, where the parts are joined. The exterior 
solutions, on the other hand, correspond to the exterior part of particular Penney solutions, which then determine 
the asymptotic properties of the solutions, such as their mass and their scalar charge. 

It appears interesting to also consider monopolc-antimonopole pairs with large electric charge. When electric charge 
is added to a monopolc-antimonopole pair, the monopole and antimonopole experience a repulsive force and the poles 
move further apart [8(. More importantly, however, the pair begins to rotate about its symmetry axis with an angular 
momentum J, equal to the product of the total electric charge Q and the (magnitude of the) individual magnetic 
charge n of the constituent magnetic poles, J = nQ [Io| . 

Preliminary study shows, that the presence of electric charge also leads to bifurcations of the charged monopole- 
antimonopolc solutions, with new branches of solutions arising. Considering larger values of the charge, we observe, 
that the values of v, cr and fh appear to be consistent with the respective Penney relations. This seems surprising 
though, since the Penney solutions carry no angular momentum, while the charged dyons possess angular momentum 
J = nQ. 
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